ASYMPTOTICS BEHAVIOUR IN ONE DIMENSIONAL MODEL OF 

INTERACTING PARTICLES 



RAFAL CELINSKI 

Abstract. Wc consider the equation ut — eu^x + {u K' * u)x ior x E M., t > and with 
e > 0, supplemented with a nonnegative, integrable initial datum. We present a class 
' of interaction kernels K' such that the large time behaviour of solutions to this initial 

CN ■ 

' value problem is described by a compactly supported self-similar profile. 

1. Introduction 

P_i I We study the asymptotic behaviour of solutions to the one-dimensional initial value 

. problem 

^ ■ (1.1) Ut = euxx + {u K' * u)^ for X G M, t > 0, 

^ (1.2) u{x, 0) = Uo{x) for x G M, 

where the interaction kernel K' is a given function, an initial datum uq G L^(]R) is 
CN ■ nonnegative and e > 0. 

00 ! Equation fll.ip arises in study of an animal aggregation as well as in some problems 

^ ! in mechanics of continuous media. The unknown function u = u{x,t) represents either 
O ■ the population density of a species or, in the case of materials applications, a particle 
density. The kernel K' in (11. II) can be understood as the derivative of a certain function 
K, that is, K' stands for dK/ dx. We use this notation to emphasise that the cell 
^ ! interaction described by equation (II. ip takes place by means of a potential K. Moreover, 
• our assumptions on interaction kernel K' imply that equation (11.11) describe particles 
interacting according to a repulsive force (this will be clarified bellow). 

Let us first notice that the one-dimensional parabolic-elliptic system of chemotaxis 

(1.3) Ut = eu^x - {uvx)x, -v^x + v = u, xgM, t>0 
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can be written as equation (11. ip . Indeed, if we put K{x) = — |e~'^' into the (11. ip . which 
is the fundamental solution of the operator d^—ld, one can rewrite the second equation 
of (II -Sp as w = —K * u. Here, however, we should emphasise that, below we consider 
repulsive phenomena, where the interaction kernel has the opposite sign, see Remark 12.31 
for more details. 

This work is motivated by the recent publication by Karch and Suzuki [S] where the 
authors study the large time asymptotics of solutions to (ll.ip -f lL2|l under the assumption 
K' G L^(M). They showed that either the fundamental solution of the heat equation or 
a nonlinear diffusion wave appear in the asymptotic expansion of solutions as t — ?■ oo. 
Analogous results on the solutions to the one dimensional chemotaxis model (II. 3p can be 
found in [121 US]- Here, we would like to point out that, in all those results, a diffusion 
phenomena play a pivotal role in the large time behaviour of solutions to problem (II. ip - 

(O. 

The main goal of this work is to show that for a large class of interaction kernels 
K' G L°°(M)\L^(M), the diffusion is completely negligible in the study of the large time 
asymptotics of solutions. Let us be more precise. Our assumption on the interaction 
kernel imply that K'(x) is sufficiently small perturbation of the function —^H(x), where, 
A G (0, oo) is a constant and H is the classical sign function given by the formula: 
H{x) = —1 for X < and H{x) = 1 for x > (c/. Remark l2.2l) . Under these assumptions, 
we show that for large values of time, a solution of problem (ll.ip -f lL2|) looks as a compactly 
supported self-similar profile, defined as the space derivative of a rarefaction wave, i.e. 
the solution of the Riemann problem for the nonviscous Burgers equation Ut + AuUr^ = 
(see Corollary 12.61 for the precise statement). 

In our reasoning, first, we consider e > 0, and our result on the large time behaviour 
are, in some sense, independent of e. Next, we pass to the limit e — )■ to obtain an 
analogous result for the inviscid aggregation equation Ut — {u K' * u)x = 0. In particular, 
our assumptions imply that weak, nonnegative solutions to the initial value problem for 
this inviscid equation exists for alH > 0. 

To conclude this introduction, we would like to recall, that the multidimensional inviscid 
aggregation equation — V ■ {uWK * n) = was derived as a macroscopic equation from 
the so-called "individual cell-based mode" [H [15], namely, as a continuum limit for a 
system of particles Xk(t) placed at the point k in time t and evolving by the system of 
differential equations: 



dt 
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where K is the potential. Results on the local and global existence as well as the blow- 
up of solutions of this inviscid aggregation equation one can find in [H El |3l [TT] and in 
references therein. 

Notation. In this work, the usual norm of the Lebesgue space L^(R) with respect to the 
spatial variable is denoted by || ■ ||p for any p G [1, oo] and 14^'^'*'(M) is the corresponding 
Sobolev space. The set C^(M) consist of smooth and compactly supported functions. 
Moreover (/ * g){x) denotes the usual convolution, i.e. (/ * g){x) = f{x — y)g{y) d|/. 
The letter C corresponds to a generic constants (always independent of x and t) which 
may vary from line to line. Sometimes, we write, e.g. C = C(q!,/5,7, ...) when we want 
to emphasise the dependence of C on parameters a, /9, 7, .... 



We begin our study of large time behaviour of solution by recalling that, for e > 0, the 
initial value problem f ll.ip - fll.2p is known to have a unique and global-in-time solution for a 
large class of initial conditions Uq and interaction kernels K'. Such results are more-or-less 
standard and the detailed reasoning can be found in [9]. In particular, our assumptions 
(see Theorem 12. II below) imply that K' G L°°(]R), hence the kernel K' is mildly singular 
in the sense stated in [9l Thm 2.5]. In this case, results from [9] can be summarised as 
follows: for every Uq G L^(R) such that Uq > 0, there exists the unique global-in-time 
solution u of problem fll.ll) - fll.2p satisfying 



In addition, the condition Uq{x) > implies u{x, t) > for all x G M and t > 0. Moreover 
we obtain the conservation of the L^-norm of nonnegative solutions: 



In Theorem 12.51 below, we pass to the limit £ -> 0, to obtain nonnegative weak solutions 
of problem fll.ip - fll.2l) with e = 0, for which the conservation of mass (12.11) holds true, as 
well. 

The goal of this work is to study the large time behaviour of solution to (ll.ip - (ll.2l) . 
First, we state conditions under which these solutions decay as t — )■ 00. 

Theorem 2.1 (Decays of norm). Assume that u = u{x,t) is a nonnegative solution 
to problem (ll.ll) - (ll.2p with e > 0, where the interaction kernel has the form K'{x) = 
— ^H{x) + V{x), where H is the sign function, A > is a constant, and the function V 



2. Main results 



uec ([0, +00), L\m)) n c ((o, +00), w^^'i(M)) n ((o, +00), L\m)) . 



(2.1) 
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satisfies 

(2.2) V e W^'^R) with IlKllii < A. 

Suppose also that uq G L^(]R) is nonnegative. Then for every p G [1, oo] the following 
inequality hold true 

(2.3) \\u{t)\[p <{A- IlKlli)"? \\uo\\y^ t'-v" 
for all t > 0. 

Remark 2.2. Notice that, under assumption (12. 2p . we have V{x) = J^^Vy{y) dy. Hence, 
we get immediately that V G L°°(]R) fl C(M), \im\x\-^ooV{x) = 0, and the following 
estimate, ||V||oo < ||K,-||i < ^ hold true. Consequently, our assumption on the interaction 
kernel K' imply that K' + G Co(]R) (continuous and decaying at infinity functions). 
This means that the kernel K' has to jump at zero exactly as the rescaled sign function 
— ^H and has to converge to the constants ±^ as x — =Foo, respectively. In some sense, 
this means that the potential K(x) looks as — at x = and as — )■ +oo. 

Remark 2.3. Our assumptions on the kernel K'{x) imply that interactions between parti- 
cles are similar as in the chemorepulsion motion in chemotaxis phenomena, namely, when 
regions of high chemical concentrations have a repulsive effect on particles. Such a model 
was studied for example in [S]. 

In the next step of this work, we derive an asymptotic profile as t — )■ oo of solutions 
(ll.ip - (ll.2l) . First, notice that if the large time behaviour of a solution to problem (11. ip - 
(II. 2p is described by the heat kernel or the nonlinear diffusion wave (as e.g. in [8]) then 
we expect the following decay rate ||M(t)||p < C t '^p for all t > 0. Observe, that the 
function u from Theorem 12.11 decays faster, hence, its asymptotic behaviour as t — )■ oo 
should be different. 

From now on, without loss of generality, we assume that u{x, t) dx = Uo{x) dx = 1. 
Indeed, due to the conservation of mass (12. ip . it suffices to replace u in equation (ll.ip by 
and K' by K' J^uodx. 
Next, let us put 

(2.4) U{x,t) = [ u{y,t)dy-\, 

where u{x,t) is the solution of (ll.ip - (ll.2p . Since u = U^, using the explicit form of the 
kernel K' (cf. Lemma [3.11 below), we obtain that the primitive U = U{x,t) satisfy the 
following equation 

(2.5) Ut = eU,, - AUU, + U,V* U,, 
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which can also be considered as a nonhnear and nonlocal perturbation of the viscous 
Burgers equation. 

Our main result says that the large time behaviour of U is described by a self-similar 
profile, given by a rarefaction wave, namely, the unique entropy solution of the Riemann 
problem for the scalar conservation law 

(2.6) W^^ + AW^W^ = 

(2.7) W''{x,0) = ^H{x). 

It is well-known (see e.g. [B]) that this rarefaction wave is given by the explicit formula 



(2.8) W^{x,t):=< 



' 1 
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Theorem 2.4 (Convergence towards rarefaction waves). Let the assumptions of Theo- 



rem \2.1\ hold true. Assume, moreover, that a nonnegative initial datum uq{x) satisfies 

(2.9) / UQ{x)dx = l, and I Uq{x)\x\ dx < oo . 

Jr Jr 

Then, there exist a constant C > independent of e such that for every t > and each 
p G (1, oo] the following estimate hold true 

(2.10) \\U{;t) - W^{;t)\\p < Ct"'^^'^'^) (log(2 + t))^(^+^) , 

where U = U{x,t) is the primitive of solution of problem fll.ip - fll.2p given by 02. 4p and 
= W^{x,t) is the rarefaction wave given by fl2l 



Next, we show that the asymptotic formula fl2.10p holds also true for weak solutions of 
problem with e = 0. 

Theorem 2.5. Assume that the kernel K' has properties stated in Theorem \2.1\ and the 
nonnegative initial condition uq G L^(]R) satisfies (12.91) . Then the initial value problem 

(2.11) Ut = -AUU^ + U.,V*U, 

(2.12) U{x, 0) = Uo{x) = I u,{y) dy-^ 
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has a weak solution U G C(R x (0, oo)) such that Ux G lv^^((0, oo), L°°(R)) that satisfies 
problem (12.111) - (12.121) in the following integral sense 

- / / Uiptdxdt- Uo{x)ip{x,0) dx = — / / U'^ip^ dx dt+ / {Vx*U)ip dx dt 
Jo Jk Jr 2 Jq J]g Jo Jr 

for all if G C^(]R x [0, +oo)). This solution satisfies 

(2.13) \\U{-,t) - W\Mp < (log(2 + t))^(^+i) , 

for a constant C > 0, for all t > 0, and each p E (1, oo]. 



Next, we use the result from Theorems 12. 4l and l2.5l to describe the large time asymptotics 
of solutions to problem (ll.ip - (ll.2p . 

Corolary 2.6. Let the assumptions either of Theorem \2.1\ or Theorem \2.5\ hold true. 
For the solution u = u{x,t) of problem (ll.ll) - (ll.2p with e > we define its rescaled 
version u^{x,t) = Xu{Xx, Xt) for A > 0, x G M and t > 0. Then, for every test function 
if G C^iM) and each to > 

uHM.)i.^-(W(.,U)^A.)i. as A^H-oo. 

Jr 

In other words, for each to > 0, the family of rescaled solutions M'*'(x,to) = Am(Ax, Ato) 
to problem (ll.ip -( lL2|) with £ > converges weakly as A — oo to the compactly supported 
self-similar profile defined as 

1 



(2.14) (W^^)^(x,to): = 



for 
At 





At 


x\ < 
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x\ > 






Y 



for 
3. Large time asymptotics 

In this section, we prove all results stated in Section [2l We begin by an elementary 
result. 

Lemma 3.1. Let H be the sign function. For all ip G ^^"'^'^(IR) the following inequality 
hold true: H * (fx = 2^9. 

Proof. First, we assume that Lp G C^(]R). Then 

P PX POO 

H*ipx= H{x-y)ipy{y)dy= (py{y) dy - ipy{y) dy = 2lp{x). 
Jr J-oo Jx 

The proof for general ip G ly^'^(M) is completed by a standard approximation argument. 

□ 



Now, we are in a position to prove Theorem 12.11 concerning the decay of solution in the 
L^-spaces. 
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Proof of Theorem \2.1[ Note, that, by (12. ip . we have ||M(t)||i = ||mo||i which imphes (12. 3p 
for p = 1. Hence, we can assume that p > 1. 

We multiply equation (II. ID by pu^'^ (recall that u is nonnegative) , integrate with 
respect to x over M, and integrate by parts to obtain 

j- [ uPdx = -^^^^^ [ [{uP/')fdx + {p-l) [ u^K'^u^dx. 

First term on the right-hand side (containing e > 0) is obviously nonpositive, hence, we 
skip it in our estimates. Using the explicit form of the kernel K' = —^H + V and Lemma 
I3.lt we rewrite the second term as follows: 

(3.1) (p-l) [ uPR' *u^dx = {p- 1) (-A [ dx+ I vFV^*u dx] . 

Jr \ Jr Jr J 

Notice, that a simple computation involving the Holder and the Young inequalities leads 

to the estimates 



(3.2) 



Vr * M dx 



< ||K * mIIp+iIIm^IIp+i < 



Hence, using (13. ip and (13. 2 p we get 

(3.3) I £ u{x, ty dx<ip- 1) {-A + 1 1 1) Mt) I nil 

Moreover, it follows from the Holder inequality (with the exponents p and ^^) that 
/ dx = / UP u p dx < ( / udxy( / u^"^^ dx j , 

which means 

(3.4) [ uP+Ux> llMollpf / uPdx)~\ 

Jr ^ Jr ^ 

because ||M(t)||i = ||ito||i- Applying estimate (13. 4p to (13.30 . we obtain the following differ- 
ential inequality for dx: 

(3.5) j^J^u{x,tydx < (p- 1) {-A + \\VM Ikolk ^ (Z''^^'^^'''^^) • 
It is easy to prove that any nonnegative solution of the differential inequality 



with a constant D > 0, satisfies 
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Hence, it follows from fl3.5p and from the assumption ||V^||i < A that 

(3.6) \\ump < (A - mu)"^ wuoWi^^ t"^ 

for all t > 0. Finally, passing to the limit p — !■ oo in (13.61) we obtain 

Mt)\\oc<{A-\\v^\\iy' 



for all t > 0. This completes the proof of Theorem 12. 1[ □ 

Let us now recall some result on smooth approximations of rarefaction waves, more 
precisely, the solution of the following Cauchy problem: 

Zt - sZ,, + AZZ, = 0, 
Z{x,0) = Zo{x) = \H{x). 

where A > 0. 

Lemma 3.2 (Hattori-Nishihara [7]). Problem (13 .71) has a unique, smooth, global-in-time 
solution Z{x,t) satisfying 

i) -1/2 < Z{x,t) < 1/2 and Z^{x,t) > for all {x,t) G M x (0,cx)); 
a) for every p E [l,oo\, there exists a constant C = C{j>) > Q independent of e > such 
that 

\\z.mp<ct-'^'/p 

and 

\\z{t)-w^{t)\\p < 



for all t > 0, where W^{x,t) is the rarefaction wave given by formula (12. Sp . 

Sketch of the proof. All results stated in Lemma ESI can be found in [7j with 
some additional improvements contained in flUl sect. 3], and they are deduced from an 
explicit formula for smooth approximation of rarefaction waves. Here however, we should 
emphasise that the authors of [7j consider equation (13. 7p with e = 1 but, by a simple 
scaling argument, we can extend those results for all e > 0. Indeed, we check that the 
function f{x,t) = Z{ex,et) satisfies ft — fxx + Aff^ = 0. Hence, by the result from [7j 
we have 

ii/.(t)iip < ctV and \\f{t) - w'm, < ct-('-y^y\ 

Now, coming back to original variables, we have 

e'^WZxi^et^pKC (ety-^e'T^ 
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and so, defining the new variable t = et, we obtain ||Z2,.(t)||p < C t p with a constant 
C independent of e. A similar reasoning should be applied in the case of the second 
inequality in Lemma I3.2IIII1 □ 
Next, we study the large time asymptotics of U{x, t) = u{y, t) dy — |, which satisfy 
equation fl2.4p . Recall that u = U^- In the proof of Theorem 12.41 we need the following 
auxiliary result. 

Lemma 3.3. Let uq satisfy conditions (12. 9p . Assume that U = U{x,t), defined by 
f l2.4p . is the solution of equation fl2.5p supplemented with the initial condition Uq{x) = 
Jloo ^o(y) dy — 1/2 and Z = Z{x,t) is the smooth approximation of the rarefaction wave, 
namely, the solution of problem f l3.7p . Then, for every tg > we have 

^^PT-7^JU{t)-Z{t)h<oo 
t>to iOg[Z + t) 

Proof. At the beginning, let us notice that assumption 02.91) on uq imply that Uo{x) G 
L^{-oo, 0) and Uo{x) - 1 G ^^(0, oo). Hence, we have that Uq ~ Zq e L^{R). 

Denoting R = U — Z and using equations fl2.5p and fl3.7p . we see that this new function 
satisfies 

We multiply this equation by sgni? (in fact, by a smooth approximation of sgni?) and 
we integrate with respect to x to obtain 



d 
di 



I \R\dx = e I i?^^sgn i? dx — ^ j (t/^ — Z^)^sgn _Rdx + j UxV *Ux^g'n.Rdx. 
Jm. Jr 2 Jig j]g 



The first term on the right-hand side of the above equation is nonpositive because this 
is the well-known Kato inequality. The second term is equal to because of the following 
calculations: 



hu'^ - Z'^)xSgnRdx = / {R^ + 2RZ)^sgnRdx 
Jr Jr ^ 



2Rx\R\dx+ / 2ZRxSgnRdx + / 2Zx\R\dx 
Jr 

-2 / ZjR\dx + 2 / ZJR\dx = 



since J^Rx\R\dx = 0. Moreover, using the Young inequality, we have 



Ux V * Ux sgn R dx 



< \\UxV*Ux\\l < ||f/J|oo||V^||l||t/x.||l. 
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Hence, by the fact that Ux{t) = u(t) and using the decay estimates from Theorem 12. II for 
p = 1 and p = oo we get the following differential inequality 



which completes the proof of Lemma I3.3[ □ 
Now, we are in a position to prove our main result about convergence the primitive of 
u towards a rarefaction wave. 



Proof of Theorem\2.4\ Let Z = Z{x,t) be the smooth approximation of the rarefaction 



wave from Lemma I3.2[ Denote R = Z — U. Hence, by Lemma 13.21 and Theorem 12. ![ we 
have 

\\RAt)\\oo = \\US) - ^x(t)||oo < h(t)||oo + ||^x.(t)||oo < C 

for a constant C > 0. Moreover, using the Sobolev-Gagliardo-Nirenberg inequality 

II Rii ^ r^ii f? II ^^'""^•^ II RiP^^"*"^^ 

1 1 1 1 P _ '-'II -^X 1 1 oo 1 1 1 1 1 5 

valid for every p G (1, C)o] and Lemma [3.31 we have 

\\u{t) - zmp < cr^'-^) (iog(2 + t))^(^+i) 

for all t > 0. 

Finally, to complete the proof, we use Lemma IX^ to replace the smooth approximation 
Z{x,t) by the rarefaction wave W^{x,t). □ 

The proof of Theorem 12.51 relies on a form of Aubin-Simon's compactness result that 
we recall below. 

Theorem 3.4 (^31 Theorem 5]). Let X , B and Y be Banach spaces satisfying X d B C. 
Y with compact embedding X G B. Assume, for 1 < p < +oo and T > 0, that 

• F is bounded in L^{0,T] X) , 

• {dtf : f eF} IS bounded m 1^(0, T;Y). 

Then F is relatively compact in L^(0,T; B) [and in C(0,T; B) if p = +oo). 

Proof of Theorem \2.5i We denote f/^ as a solution of equation fl2.5p with e > supple- 
mented with a initial condition (12.121) . The proof follows three steps: first we show that 
the family 

J^ = {U' -.ee {0,1]}, 

is relative compact in C{[ti,t2],C[—R, R\) for every < ti < t2 < oo and every R > 0. 
Next, we show that there exist a function U = limg^o U'^ which is a weak solution of 
problem fl2.lip - fl2.12p . Finally we prove that U satisfy estimate fl2.13p . 
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Step 1. Compactness. We apply Theorem 13.41 with p = oo, F = and 



X = C\[-R,R]), B = C{[-R,R]), Y = W~'^\[-R,R]), 



where i? > is fixed and arbitrary, and Y is the dual space oiW^'\[-R,R]). Obviously, 
the embedding X C i? is compact by the Arzela-Ascoli theorem. 

First, we show that the sets J-" and {dxU'^ : e G (0,1]} are bounded subsets of 
L°° {\ti^t'2\.,C{[—R^ R\)). Indeed, it follows from definition of function U^, namely from 
(12:41). that 



\U'{x,t)\ < \\{U%{;t)\\i + 



\uo\\i + 



1 



Moreover, using Theorem 12.11 we have 



(3.9) 



\\{u%{-Moo<{A-m\,)-u 



1 



To check the second condition of Aubin-Simon's compactness criterion, it is suffices 
to show that there is a positive constant C which independent of e G (0, 1] such that 
sup^gj^j^^^] ll^^tf^'^lly ^ C. Let us show this estimate by a duality argument. For every 

G {{-R,R)) and t G [^1,^2], by (jSSD, (jM]) and Theorem EH we have 



dtU%t)^ dx 



< 



< \Wx\ 



\um^^+A\\u'moo\Mc 



AU'{t)Ul{t)^dx 
. I \Ul{t)\ dx+\\U, 



umv*ui{t)^dx 



m\\io\ 



v\\iM 



< 



|mo||i + A||mo||i(||mo||i + l/2)|b||oo + {A- IlKlli) ti ll^llilblli. 



Hence, the proof of Step 1 is completed. 

Step 2. Limit function. By Step 1, for every < ti < ^2 < +00, the family {If^ : e G 
(0,1]} is relatively compact in C{[ti,t2],C{—R, R)). Consequently, by a diagonal argu- 
ment, there exists a sequence of {f/^" : En G (0, 1]} and a function f7 G C((0, +00), C(M)) 
such that 



(3.10) 



U'" -> f/ as En^O in L^^(R x (0, +00)). 



Moreover, by the Banach-Alaoglu Theorem, it follows from the estimate (13.91) that 



Ux" Ux as En^O 



weak-* in L~ ((0, 00), L°°(R)) . 
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Now, multiplying equation f l2.5p by a test function G C^(Mx [0, +C)o)) and integrating 
the resulting equation over M x [0, oo), we obtain the identity 

(3.11) - / / f/""(^tdxdt- / f/o(x)(^(x,0)dx = e„ / / dx dt 

Jo Jr Jr Jo Jr 

f oo P POO 



A 
+ 2 







poo p 

(f/"") dxdt+ / f/f (K- * f/'") dx dt 
Jo Jr 



It is easy to pass to the limit — > in left-hand side of (13. lip , using the Lebesgue 
dominated convergence theorem. To deal with term in the right-hand side we make the 
following decomposition: 

(3.12) [ U^- (K *U'")^dx= I f/^" * {W- -U))^dx+ I Ul^ *U)^ dx. 
Jr Jr Jr 

We can estimate the first term on the right-hand side of (13.120 as follows: 



(3.13) 



/ f/^- (K. * [U'" - U))y^dx < m-{t)\U [ |K. * (f/'" - U)v\ dx 
Jr Jr 



Let us notice, that * {If^" — U) tends to zero as e„ — t- by Lebesgue dominated 
convergence theorem and it is bounded independently of £„. Hence, using the Lebesgue 
dominated convergence theorem and Theorem 12.11 we deduce that the right-hand side 
of (I3.13P converge to zero. The second term on the right-hand side of (I3.12p obviously 
converge to j^U^ {Vx * f7) v^dx by the weakly-* convergence of U^" in L°°(M) since {Vx * 
U)ip G L^(]R). This completes the proof of Step 2. 

Step 3. Convergence towards rarefaction wave. To prove (I2.13p . we use the Fatou 
Lemma and (I3.10p . to obtain 

\\m - W^m, < hminf ||[/-(t) - W\t)\\, 

for all t > 0. 

Now, it is enough to use Theorem 12.41 to estimate the quantity on right-hand side, since 
constant C in (I2.10p is independent of e. Hence the proof of Theorem 12.51 is finished. □ 



At last, we prove Corollary [2? 

Proof of Corollary \2.b\ First, we express the result stated in Theorems 12.41 and 12.51 in 
another way. We consider the rescaled family of function U^{x.,t) = U{Xx,Xt) for all 
A > 0. Let us also notice that W^{x, t) is self-similar in the sense that {W^) {x,t) = 
W^{x, t) for all a; G M, t > 0, A > 0. Hence, changing the variables and using Theorem 12.41 
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and Theorem 12.51 for the case e = 0, we obtain 

\\U\;to) - (W^)' (-,^0)11^ = A-i/^||f/(-, Ato) - W\;Xto)\\, < 

CA-^/P(Ato)"^^^"^) (log(2 + Ato))^^^^^) ^ 

as A — 7- oo. It means that the family of functions converge in ^^(M) as A — t- oo towards 
W^{x, t) for every to > and p G (1, C)o]. 

This scahng argument allows us to express the convergence of solutions to original 
problem f lLip -f lT^ towards a self-similar profile. Indeed, let us note that since u = U^, it 
follows immediately that u^{x, t) = Xu{Xx, At) = dxU'^{x, t). Hence, the weak convergence 
of u'^ towards the distributional derivative of the rarefaction wave d^W^ is the immediate 
consequence of the Lebesgue dominated convergence theorem and of Theorem 12.41 for 
p = oo since |?7'^(x, to)| < I^'^oix) dx + ^ □ 
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